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Abstract 

An anomaly in the distribution of quasar magnitudes based on the SDSS survey, has 
been recently reported by Longo (2012). The angular size of this anomaly is of the order of 
±15° on the sky. A low surface brightness smooth structure in 7-rays, coincides with the 
sky location and extent of the quasar anomaly, and is close to the Northern component of a 
pair of 7-ray bubbles discovered in the Fermi Gamma-ray Space Telescope survey. Molecular 
clouds are thought to be illuminated by cosmic rays. A large fraction of molecular gas in the 
Galaxy is under the form of cold H2, and may be a significant component of dark matter. 
I test the hypothesis that the magnitude anomaly in the quasar distribution, is due to a 
lensing effect by an hypothetic Supergiant Molecular Cloud (SGMC) in the Galactic halo. 
A series of grid lens models are built by assuming firstly that a SGMC is a lattice with 
clumps of 10~ 3 M Q , 10 AU in size, and assuming various filling factors of the cloud, and 
secondly a fractal structure. Local amplifications are calculated for these lenses by using 
the public software LensTool, and the single plane approximation. A complex network of 
caustics due to the clumpy structure is present. Our best single plane lens model capable of 
explaining Longo's effect, at least in sparse regions, requires a mass (1.5 — 4.1) x 10 10 M 
within 8.7 x 8.7 x (5 — 8.6) kpc 3 at a lens plane distance of 20 kpc. It is constructed from 
a molecular cloud building block of 5 x 10 5 M Q within a scale of 30 pc expanded by fractal 
scaling with dimension D = 1.8 — 2 up to 5-8.6 kpc for the SGMC. If such a Supergiant 
Molecular Cloud were demonstrated, it might be part of a lens explanation for the luminous 
anomaly discovered in quasars and in red galaxies. The mass budget may be varied by 
changing the cloud depth and the fractal dimension. 



1 Introduction 



An anomaly in the distribution of quasar magnitudes based on the SDSS survey, was recently 
reported by Longo (2012). The effect, which is an enhancement of the order of 0.2 mag, is 
statistically extremely significant in amplitude, and rather well defined in angular size. A 
first hint of this anomaly was in fact visible in quasar data as early as three decades ago 
(Giraud & Vigier 1983). Abate & Feldman (2012) reported a similar enhancement in the 
distribution of SDSS luminous red galaxies. Because the luminous anomaly does not seem 
to depend on color, a possible physical explanation is gravitational lensing. Nevertheless, 
if the lens is at cosmological distance, its huge angular extent translates into a lens radius 
~ 350 Mpc, and its critical surface density implies a total mass of the order of 10 21 M Q , which 
is barely consistent with ACDM cosmology. A major reason to test a lensing hypothesis with 
a lens at short distance from the observer, possibly in the Galactic halo, is the large angular 
extent of the effect combined with its low amplitude. 

Two large 7-ray low surface brightness bubbles, extending 50° above and below the 
Galactic center were discovered by Su et al (2010) in Fermi Gamma-ray Space Telescope 
survey. These bubbles are spatially correlated with a microwave excess in WMAP. Both the 
Fermi bubbles and the WMAP haze can be explained by synchrotron and inverse- Compt on 
from a hard electron population (Dobler et al 2010). In addition to the bubbles in the Fermi 
map, there is a single large faint structure East of the Northern bubble (I ~ —60°, b rs 60°), 
that is at the location of the quasar anomaly. Any molecular cloud, if present in this faint 
structure, would be illuminated by cosmic rays escaping from the Northern bubble. Cosmic 
ray/nucleon interactions in molecular clouds were discussed by e.q. Aharonian & Atoyan 
(1996), de Paolis et al (1999), Kalberla et al (1999), Sciama (1999), Shchekinov (2000). The 
spectra of 7-rays arising from molecular clouds illuminated by cosmic rays were recently re- 
examined by Ohira et al (2011). The angular extension of the 7-ray structure is consistent 
with that of the quasar luminosity anomaly. We make the working hypothesis that 7-rays 
coming from the faint structure trace a SuperGiant cold Molecular Cloud (SGMC) in the 
Galactic halo. 

A large fraction of molecular gas in the Galaxy is under the form of cold H2, and may 
be a significant component of dark matter (Pfenniger, Combes, & Martinet 1994). The 
variation of universal rotation curves (Persic, Salucci, & Stel 1996) with mass in spiral 
galaxies may be interpreted as an increasing fraction of baryonic dark matter with decreasing 
mass of galaxies (Giraud 2000a,b; 2001). The fractal structure of molecular gas clouds, their 
formation, thermodynamics and stability were intensively studied (Scalo 1985; Henriksen 
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1991; Falgarone, Phillips, & Walker 1991; Walker & Wardle 1998; Irwin, Widraw, & English 
2000). The existence of small dark clumps of gas in the Galaxy was invoked to explain 
Extreme Scattering Events of quasars (Fielder et al 1987, 1994). The duration of such 
events suggests sizes and masses of individual cloudlets of the order 10 AU and 10~ 3 M . 
These values of mass and size are also allowed by the searches for lensing toward the LMC 
(Draine, 1998). Cloudlets can survive if the number of collisions is sufficiently small and their 
stability is best explained if they are self-gravitating (Combes 2000). Concerning lensing, 
the average surface density of a SGMC would be sub-critical, but the chaotic structure of 
the matter distribution may lead to a network of caustics with high local magnification. 

In the present Letter I test the hypothesis that the magnitude anomaly in the quasar 
distribution can be explained by the lensing effect of an hypothetic SGMC in the Galactic 
halo. 

2 Micro lensing and cloud lens models 
2.1 Introduction 

Lensing of distant galaxies and quasars shears and magnifies their images. There are ex- 
tensive reviews and books on lensing including Schneider, Ehlers & Falco (1992), Mellier 
(1999), Schneider (2005, extensive references therein) and public algorithms. Fundamental 
theory, ray-tracing algorithms (Valle & White 2003 and references therein), and observational 
techniques have been developed for two to three decades to infer cluster mass distributions 
(review by Gladders et al 2002), large scale structures and cosmological models (Soucail et 
al 2004; Wambsganss et al 1997, 2004; Sand et al 2005). 

Systematic studies of micro- lensing were initiated in the early 90s, by Wambsganss (1990), 
who developed ray-shooting simulations to show the magnification pattern of star fields. In 
the present paper we will mainly use the public software LensTool (Julio et al 2007) which 
includes strong and weak lens regimes. 

The change in direction of a light ray can be written 

da = —2 grad (f> d\ 

where da is the bend angle, grad is the spatial gradient perpendicular to the light path, 
is gravitational potential, and x is the radial co-moving coordinate. The distortion can be 
entirely described by the linearized lens mapping through the Jacobi matrix A of transform. 
This matrix is normally expressed with three terms depending on derivatives of the the 
lensing potential: the convergence k, the amplitude of the shear 7 and the rotation u>. The 
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convergence changes the radius of a circular source, the shear changes its ellipticity, the 
rotation angle is the phase of the shear. Rotation occurs in the case of multiple lens planes. 

The convergence k = \{d\ + at location 9 is related to the projected surface density 
E(0) at that location by 

^crit 

in which the critical surface density E crit is given by 

c 2 D a 
crit ~ AixG DiDi s 

where D s , D% and D ts are respectively the distances observer to source, observer to lens, and 
lens to source. 

The components of the shear are given by 71 = \{d\ — <9f ) and 72 = d\di 0. The shear 
amplitude is the module of its two components. 

The present Letter focusses on micro-lensing by numerous clumps; the convergence and 
the shear of any underlying smooth component will be ignored. 

2.2 Preliminary test lens models 

The critical surface density of a smooth lens located in the Galactic halo on distant objects 
such that D s ps Di s and assuming Di = 15 kpc is 10 5 M pc~ 2 , which leads to a mass of 
~ 10 14 M for a kpc scale SGMC, so a 3D filling factor of the order of 1CT 5 is required to 
fit in the halo. A series of cloud grid models of kpc size in depth, with individual clumps of 
mass 1CT 3 M , 10 AU in diameter, are built with various filling factors. 

Cloud models are constructed as follows: an initial fits frame of 4000 x 4000 pixels with 
zero intensity including grids of the order of 10-100 pixels with common value 1 is created, 
so that typical separations of the order of A ~ 10000 AU are reached by assuming a physical 
pixel size of 10 AU. A series of layers is obtained by shifting randomly the coordinates of the 
initial frame which we co-add. We repeat this shift and add procedure on the first co-added 
layer to get a second generation, and repeat again the process up to several generations. For 
example a model with surface density 0.8 E crit and 1 kpc in depth, is built with 8 5 + 8 4 layers, 
from an initial frame with 55 non-zero pixels. This lens has 2 x 10 6 micro-lenses, some being 
combined on the same pixel. 

Lens mass models are obtained by multiplying the final frames by the mass of an individ- 
ual clump, and transformed into ascii files assuming that luminous pixels are point masses. 
Point masses are used for simplicity. A detailed lens analysis, in which self-gravitating gas 
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clumps are described by polytropic models, can be found in Draine (1998). Lens simulations 
with LensTool are done on small sections of the final frames including a sufficiently large 
number of point masses each. 
Three cases are considered: 

- ac — 0.2 in order to compare our simulations with those in the literature, 

- Ac — 0.09 because the theoretical mean amplification calculated from 



and assuming 7 = is that of Longo's effect < [i >= 1.2, 

- k = 0.8 and 0.9 to explore cases with high magnification and dense caustics. Parameters 
of test lens models are given in Table 1. 

In the case k = 0.2, individual caustics can be identified, there are perturbations of single 
point caustics by other point masses giving complicated caustic structures comparable with 
those shown in Wambganss (1990). The average amplification value obtained by simulations 
is the same as the analytic value. 

For k = 0.09, the magnification pattern differs from the case k = 0.2: individual caustics 
can be identified but there is almost no perturbation of single clump caustics by other 
clumps presumably because the density is low. So there is no complicated caustic structure 
comparable to the previous case. The average amplification value is the same as the analytic 
value with a dispersion due to the small number of clumps per frame section. 

For k = 0.8, the magnification pattern shows a high clustering of caustics. While the 
average magnification in simulations, < fi >= 25, is the same as the analytic value, large 
variations are detected with some locations showing huge amplification. We recover here the 
non-linear spatial distribution of amplifications mentioned by Wambganss (1990). 

All these kpc 3 test models are by far too massive to be used as elementary modules of a 
SuperGiant Cloud of the order of 400 kpc 3 . Low filling factors are necessary. 

A model with local surface density ac = 0.8 and a 2D filling factor of / = 1.2/25 = 0.048, 
namely a series of lenses with high amplification 25 diluted in a surface ~ 20 larger, was 
considered (test model E). This patchy model gives the expected average magnification 
< n >= 1.2 with a mass 2.3 times smaller than the case k = 0.09 with / = 1. The case 
k = 0.9 and a 2D filling factor of / = 1.2/100 = 0.012 provides very high local magnification 
up to 100, an average < \i >= 1.2 for a total mass reduced by a factor 9 (test model F). 
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Table 1: Main parameters for lens mass models. The columns indicate respectively: model 
identification, mean separation A between clumps, 2 D filling factor /, average surface 
density in filled area E for a 1 kpc depth, average amplification < /i >, mass density p per 
pc 3 , mass Mi of a SGMC of 1 kpc 3 . In models E & F the density is given for the dense 
regions, and the magnification is given both for the dense regions and on average. 

Model A (AU) / g (M pc' 2 ) < p > p (M pc~ 3 ) M 1 (M Q ) 



A 


10000 


1 


8 x 10 3 


1.18 


8 


8 x 10 9 


B 


9856 


1 


9 x 10 3 


1.21 


9 


9 x 10 9 


C 


7590 


1 


2 x 10 4 


1.53 


20 


2 x 10 10 


D 


4760 


1 


8 x 10 4 


25 


80 


8 x 10 10 


E 


4760 


0.048 


8 x 10 4 


25, 1.2 


80 


3.8 x 10 9 


F 


4575 


0.012 


9 x 10 4 


100, 1.2 


90 


10 9 



2.3 Constrained lens models 

Our lens mass models must have densities lower than the limit of molecular cloud cores: if 
we assume a maximum mass ~ 10 M Q within a scale length of L ~ 0.1 pc above which 
molecular clouds would collapse (Kayama et al 1996), and extrapolate this density through 
Larson's (1981) scaling law p oc L~ L1 , one gets an upper limit density ~ 350 M pc~ 3 
for a scale length of ~ 100 pc. Whereas our densities satisfy this criterium, a molecular 
cloud in equilibrium with scale 30 pc, close to our building block, would have a typical mass 
~ 5 x 10 5 M (Solomon et al 1987), and a density 18.5 M pc" 3 which is significantly below 
that of models D, E, & F. 

Reducing the size of the grid cells quickly boosts the amplification. Therefore a solution 
may come from a fractal structure of Molecular Clouds with rather small cells satisfying 
the equilibrium condition, separated by less dense regions, all assembled in Giant Molecular 
Clouds up to a SGMC, as it was proposed by Pfenniger & Combes (1994; hereafter PC) for 
cold dark H 2 clouds. A cloud cell of 5 x 10 5 M for a scale 30 pc will be used as building 
block for fractal SGMC models below. 

If masses obey a fractal scaling relation with dimension D from mass M min to M max and 
scale L min to L max one may write: 

M f L \ D 

-^min V L mm J 

up to -L max - The mass at level I consists in a number of clumps N at level / — 1 related to 
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Table 2: Main parameters for fractal lens mass models. The columns indicate respectively: 
model identification, lens plane distance D h fractal dimension D, mass of building block 
M min , scale of building block L min in pc, scale of SGMC L in kpc, number N of building 
blocks within 125 kpc 3 , maximum integrated surface density £ , median magnification fi, 
mass M 5 of a SGMC within 125 kpc 3 . 



Model D l D M min (M ) L min (pc) L(kpc) N g (M pc' 2 ) fj, M 5 (M ) 



G 


16 


1.64 


5 x 10 5 


30 


5 


4 500 


8.9 x 10 3 


1.03 


2.2 x 10 9 


H 


16 


2 


5 x 10 5 


30 


5 


27 750 


1.7 x 10 4 


1.12 


1.4 x 10 10 


I 


20 


1.8 


5 x 10 5 


30 


5 


10 106 


1.2 x 10 4 


1.12 


5 x 10 9 



the scale by the relation: 

hhl = n~ 1/d 
Li 

Clouds can survive if the number of collisions is sufficiently small, which according to PC 
and Irwin, Widraw, & English (2000) requires 1 < D < 2. If D > 2.33, internal and external 
collisions are too numerous for maintaining a fragmented cloud. Two cloud models with 
two fractal dimensions are considered first: the case D = 1.64 which was suggested in the 
original PC paper, and the upper value D = 2. 

Extrapolating from L min = 30 pc, M min = 5 x 10 5 M Q to L max = 5 kpc, and assuming a 
fractal dimension D = 1.64, one gets a total mass M max = 2.2 x 10 9 M within 125 kpc 3 , and 
iV ~ 4500 cloud cells, i.e. ~ 16 building blocks on each axis. The surface density integrated 
over a column L max = 5 kpc of 16 blocks is E = 8.9 x 10 3 M pc~ 2 , which is the maximum. 
Parameters of this model are given in Table 2 (model G). The median magnification derived 
from the probability p(n) of a light ray to cross < n < 16 cloud blocks, is found to be 
J2 n{n)p(n) = 1.03, so half the surface provides a magnification between 1.20 and 1.03 which 
is not sufficient to explain Longo's effect. 

The case D = 2 has a different behavior (model H). Being more populated than for 
D = 1.67, the densities are higher, the median amplification is /i — 1.12, and 25 % of the 
surface has 1.18 < /j, < 1.48. Regions with ji > 1.18 do show networks of caustics and non 
gaussian amplification distributions. Model H appears at the limit of explaining Longo's 
effect with a mass of 4.1 x 10 10 M within 28°. 

2.4 Shear 

The shear was not considered in our previous models. This is due to the fact that the shear 
is a term resulting from the mass outside the light beam, like a smooth galactic component 
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to a star field. In our case however, cloud blocks are mass concentrations which may induce 
a shear component. A series of tests were done on sections of our grids using LensTool and 
were repeated with the MRLENS package (Starck et al 2006) which is dedicated to weak 
lensing. There is considerable dispersion between grid sections with the same average surface 
density. Shear is not detected on grid sections with a single block (5 x 10 5 M ), that is for 
grids with measured amplification \i = 1.01. Shear is significant in grid sections of average 
amplification ji = 1.1 and above, but there is an important spread 0.1 k < 7 < k. Assuming 
7 = k, the amplification would slightly change, so the same amplifications as for model 
with D = 2 would be obtained with D = 1.93. In that case the mass within 125 kpc 3 is 
9.8 x 10 9 M , and that of a cloud of 5 kpc in depth extending over 28° is 2.9 x 10 10 M . 

2.5 Distance to the lens plane 

Locating the lens plane at an observer distance D\ = 20 kpc rather than 16 kpc, reduces the 
critical surface density by 25% to E cr it = 7.5 x 10 4 M pc~ 2 for infinitely distant objects. 
For this critical surface density the same amplifications as in model H, are obtained for a 
cloud with D = 1.80, a mass 5 x 10 9 M within 125 kpc 3 , and a total mass 1.5 x 10 10 M 
for a SGMC of 5 kpc in depth and extending over 8.7 x 8.7 kpc 2 (model I). 

Increasing the depth to 8.6 kpc and using a dimension D = 2, we obtain model J which 
provides a median amplication of fi = 1.15, with 14 % of the surface amplihed by more 
than 1.5. The total mass total within 8.6 x 8.6 x 8.6 kpc 3 is 4.1 x 10 10 M . 

A single plane lens model has been assumed. If the cloud were spread over a large distance 
in depth, extending for example from D\ = 10 kpc to 30 kpc, the critical surface density 
would vary by a factor 3 with D h the distant part of the lens being more efficient than the 
nearby one. Over that distance, the cosmic-ray flux coming from the central part of the 
Galaxy would decrease by a factor 9, so if the cloud had a large depth, there would be a 
distance shift between the illuminated area and the average lens plane. The mass budget 
may be further optimized by changing the depth of the SGMC and its fractal dimension, 
and consequently the probability of a light-ray to cross a sufficiently large number of cloud 
blocks. 

2.6 H2 clumps and WIMPS 

H2 clumps are extremely dense compared with WIMPs clumps of the same mass, as predicted 
by N-body simulations (Diemand et al 2005, Springel et al 2008). Let consider for example 
a WIMP clump seeded by an H 2 clump of the same mass M = 10~ 3 M . The WIMP clump 
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alone would have a radius of 1 000 AU in the absence of H 2 . Under the attraction of the 
central H 2 clump, it will experience a first infall during which its density will increase by 
10 6 in ~ 1 Myr. If WIMPS are low mass neutralinos of m x ~ 10 GeV, number density 
n x , and cross-section av=o = 10~ 26 cm 3 s _1 for annihilation (Cerdeno et al 2012), those 
concentrated in the clump will have significant annihilation, leading to a mass transform 
1.4 10~ 7 M Q Myr -1 mainly into p, p— and 7. The cumulated mass loss of those neutralinos 
is 50% in ~ 6 Gyr, thus an upper limit of ~ 50 % the clump mass may remain in WIMPs 
within 10 AU over cosmological ages. 

3 Conclusion 

A luminosity excess of 0.2 mag was recently discovered in the SSDS survey of quasars. This 
effect which has a huge spatial extent of the order of 30°, is extremely well determined 
statistically, and coincides with a large and diffuse luminosity in 7 detected by the Fermi 
Gamma-ray Space Telescope. 

We make the working hypothesis that the detected 7 luminosity traces a SuperGiant 
Molecular Cloud illuminated by cosmic rays, at least partly, and test the hypothesis on 
whether such an cloud would be capable of explaining the quasar luminosity excess by 
micro-lensing. 

Cloud models are built by using compact clumps of mass 10~ 3 M Q , 10 AU in size, as 
suggested by Extreme Scattering Events on quasars, typical molecular cloud building blocks 
of mass 5 x 10 5 M for a scale of 30 pc, and a fractal structure of the SuperGiant Cloud 
with dimension 1.64 < D < 2 up to a scale of 5-8.6 kpc, and located at an observer-lens 
distance D t = 16 — 20 kpc. 

A SGMC located at a distance Di = 20 kpc from the observer with fractal dimension 
D = 1.8, total mass 1.5 x 10 10 M , 5 kpc in depth and extending over 8.7 kpc do show 
amplifications of the order ofLongo's effect in sparse regions. A SGMC with depth increased 
to 8.6 kpc and dimension D = 2 appears capable of explaining Longo's effect with a huge 
mass of 4.1 x 10 10 M . 
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